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We study quantum transport and scattering of massless Dirac fermions by spatially localized static 
magnetic fields. The employed model describes in a unified manner the effects of orbital magnetic 
fields, Zeeman and exchange fields in topological insulators, and the pseudo-magnetic fields caused 
by strain or defects in monolayer graphene. The general scattering theory is formulated, and for 
radially symmetric fields, the scattering amplitude and the total and transport cross sections are 
expressed in terms of phase shifts. As applications, we study ring-shaped magnetic fields (including 
the Aharanov-Bohm geometry) and scattering by magnetic dipoles. 
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I. INTRODUCTION 

The recent theoretical prediction and subsequential ex- 
perimental verification of the conducting surface state ex- 
isting in a strong topological insulator (TI) has generated 
a burst of activity, reviewed in Refs. [l| and y. In a TI, 
strong spin-orbit couplings and band inversion conspire 
to produce a unique time-reversal invariant topological 
state different from a conventional band insulator. Using 
Bi2Se3 as reference TI material, one finds a rather large 
bulk gap w 0.3 meV, and surface probe experiments have 
provided strong evidence for the topologically protected 
gapless surface stated The measured spin texture of the 
surface is well described by two-dimensional (2D) mass- 
less Dirac fermions, where the spinor wavefunction has 
precisely two entries corresponding to physical spin. Un- 
der this "relativistic" description, spin and momentum 
are always perpendicular, and the surface state is sta- 
ble against the effects of weak disorder and weak inter- 
actions due to the underlying topological protection.^ 
Useful insights can then already be obtained from a non- 
interacting disorder-free description. Massless 2D Dirac 
fermions are also realized in single carbon monolayers of 
graphene, for reviews see Refs. U, [5| and |g. The limit of 
ballistic transport in this 2D material seems experimen- 
tally within reachj^ and graphene experiments have re- 
ported characteristic Dirac fermion signatures, e.g., Klein 
tunneling ^^ In most experiments performed so far, cor- 
relation effects turned out to be weak, and it is again of 
interest to examine the single-particle theory in the ab- 
sence of disorder. In contrast to the TI surface case, the 
two entries of the spinor wavefunction in graphene corre- 
spond to the two atoms in the basis of graphene's honey- 
comb lattice. Furthermore, there are four Dirac fermion 
"flavors" in graphene, due to the physical spin and the 
KK 1 ( "valley" ) orbital degeneracy* 5 - Another condensed- 
matter realization of Dirac fermions is given by the quasi- 
particles in <i-wave superconductors, 11 but to be specific, 
we here focus on the surface state in a TI and on mono- 
layer graphene. 

The fact that the electronic properties of both these 
applications correspond to massless 2D Dirac fermions 
calls for a unified description of their transport proper- 



ties. In the graphene context, much theoretical effort has 
been devoted to advancing the scattering theory of Dirac 
fermions in electrostatic potentials* 5 -! 1 ^ in particular for 
the Coulomb impurity J^ In this paper, we instead study 
the scattering of massless Dirac fermions by a local mag- 
netostatic perturbation. The model, see Eq. (j2.3[) below, 
describes, in a unified manner, the effects of spatially in- 
homogeneous orbital magnetic fields, exchange-mediated 
fields due to adjacent ferromagnetic (FM) layers and Zee- 
man fields in topological insulators, as well as strain- or 
defect-induced pseudo-magnetic fields in graphene. For 
the Schrodinger fermions realized in 2D semiconductor 
electron gases, such perturbations, e.g., magnetically de- 
fined barriers, steps, and quantum wells, have been inves- 
tigated both theoretically-^— and experimentally* 1 ^ 1 ^ 
The desired magnetic field profiles were generated by de- 
position of lithographically patterned FM layers on top 
of the sample. For soft FM materials, one can change the 
magnetization orientation by weak magnetic fields. An- 
other possibility is to use a type-II superconductor film 
instead of the FM layer ^21 

Previous theory work on TIs in inhomogeneous mag- 
netic fields has addressed only a few setups. For the 
transmission of an electron through a magnetic barrier 
(assumed homogeneous in the transverse direction), as 
a function of either exchange field or applied bias volt- 
age, Mondal et al*^ predict an oscillatory behavior or 
even a complete suppression of the transmission proba- 
bility, and hence of the conductance. A spin valve geom- 
etry with two adjacent magnetic barriers, characterized 
by non-collincar exchange fields, has also been studied. 23 
As a model for a classical magnetic impurity, the spin- 
resolved density of states was calculated for a disc-shaped 
magnetic field profile ! 24 ' 25 

For graphene, a vector potential perturbation can 
again be due to external orbital fields, but may also de- 
scribe the effects of strain 5 -! 2 ^— and dislocations or other 
topological defects* 3 ^ Several theoretical works have ad- 
dressed aspects of the electronic structure and the trans- 
mission properties for Dirac fermions in graphene in the 
presence of inhomogeneous magnetic fields. The sim- 
plest case is encountered for effectively ID problems with 
translational invariance in the (say) y-direction, e.g., for 



a magnetic step or a magnetic barrier £2rM. For suit- 
able ID magnetic field profiles, it is possible to have 
magnetic waveguides (along the ^-direction) , 35 ' 36 where 
electron-electron interaction effects play an important 
roleM- Periodic magnetic fields, i.e., ID magnetic su- 
perlattices, have also been addressed.— ~— For radially 
symmetric fields, total angular momentum conservation 
again simplifies the problem and gives an effective ID 
theory. This has allowed for studies of quantum dot 
or antidot geometries , 32 i 33 ] 42 i 43 where true bound states, 
not affected by Klein tunneling, may exist. In quan- 
tum dot setups, interaction effects become important for 
strong confinement. 44 When the vector potential corre- 
sponds to an infinitely thin solenoid, we encounter an 
ultra-relativistic Dirac fermion generalization of the cele- 
brated Aharonov-Bohm (AB) calculation ^ 46 This gen- 
eralization was discussed before,— ~— and exact results 
for the transmission amplitude can be deduced. Recent 
studies have also addressed the current induced by an 
AB flux— and the behavior of the conductance when 
the chemical potential is precisely at the Dirac neutrality 
point^ Such an AB conductance can be probed experi- 
mentally in ring-shaped graphene devicesi 55 ' 56 

In this article, we formulate a general scattering theory 
approach for massless 2D Dirac fermions in the presence 
of such magnetic perturbations. In Sec. m we introduce 
the model and outline its application to graphene and 
topological insulators. In Sec. IIII1 we formulate the gen- 
eral scattering theory, previously given for electrostatic 
potentials^ for the magnetic case. The scattering am- 
plitude and cross section are specified, and we discuss the 
Born approximation in Sec. IIII Al For the radially sym- 
metric case, total angular momentum conservation allows 
to express the scattering amplitude in terms of phase 
shifts for given total angular momentum, see Sec. IIIIBI 
In sections IIVI and [Vj we present applications of this 
formalism. First, in Sec. IIV| we discuss scattering by 
a magnetic dipole within the Born approximation. Sec- 
ond, in Sec. El we consider ring-shaped field profiles. This 
case also contains the AB solenoid in a certain limit, see 
Sec. IV B( and we discuss how our phase-shift analysis re- 
covers known results for the AB effect. We also address 
the scattering resonances appearing due to quasi-bound 
states in such a ring-shaped magnetic confinement. Fi- 
nally, some concluding remarks can be found in Sec. I VII 



II. 



DIRAC FERMIONS IN GRAPHENE AND 
TOPOLOGICAL INSULATORS 

A. Model 



strain-induced pseudo-magnetic fields. In fact, for those 
applications we have A z = 0, see below. In addition, for 
the TI case, we allow for a Zeeman field or for exchange 
fields caused by nearby ferromagnets, whose components 
are contained in the field M(r) — (M x ,M y ,M z ), where 
prefactors such as the Bohr magneton or the Lande factor 
are included. With the Pauli matrices <x = (cr x ,ay,a z ) 
and the momentum operator p — —ih(d x ,d v ,0), the 
single-particle model reads (e > 0) 



H = v F cr ■ Ip + -A{r)\ + cr ■ M{r) - eV(r). (2.1) 

The Fermi velocity in graphene is vp — 10 6 m/s, while 
for a TI surface state, a typical value 3 for Bi 2 Se3 is vp « 
5 x 10 5 m/s. The low-energy Hamiltonian in Eq. (|2.1j) is 
valid on energy scales close to the neutrality level (Dirac 
point), well below the bulk band gap for the TI case and 
within a window of size ~ 0.5 eV for graphene. Both the 
vector potential A and the field M can now be combined 
to a vector field 



A(r) = A + M, 

evp 



(2.2) 



which contains all considered "magnetic" perturbations. 
Interesting physics also follows in the presence of both 
A(r) and a scalar potential V(r), but we put V(r) = 
below. 

For the formulation of the scattering theory, it is conve- 
nient to employ cylindrical coordinates, x = rcos</> and 
y = rsin0, with unit vectors e r = (cos</>, sin</>, 0), e^ = 
(— sin </>, coscj), 0), and e z . With A = A r e r + A^e^-^ A z e Zl 
Eq. (|2.1|) takes the compact form 



H = VFe-^^He 14 "** 12 , 
H = I —ihd r + -A r ) a x 4 



(2.3) 



1 T e A I 

~Jz + -Ad, (Ty, 

r c ' " 



where the total angular momentum operator is 

J z = —ihdff, + ha z /2. 



(2.4) 



For the case of azimuthal symmetry, d^A^.r — 0, this is 
a conserved quantity, [J Zl H] =0, with eigenvalues Hj for 
half-integer j. In Eq. (|2.3|) we have put A z = 0, which is 
the case for all fields studied below. 

Let us now discuss how Eq. (|2.3p relates to the sur- 
face states in a strong topological insulator, where the 
spin direction is tangential to the surface and perpendic- 
ular to momentum. For a given two-component spinor 
wavefunction ip( r )> the spin (s) and particle current (J) 
density operators are 1 ^ 



In this section, we describe the 2D Dirac fermion model 
for electronic transport in graphene or the TI surface 
studied in this work. As perturbations, we first allow for 
an external static vector potential, A(r) = (A x ,A y ,A z ) 
with r = (x,y), which is included by minimal coupling 
and describes orbital magnetic fields and (for graphene) 



s(r) = -V> f (e 2 x er)V>, 

j(r) = vpip^ (a x e x + a y e y )il;, 



(2.5) 



i.e., both spin and current are confined to the surface 
and obey s ■ j = 0. Under stationary conditions, the 



continuity equation, dt(ip ip) 
the relation 



E<=*,t, d di = °: implies 



Y^ 9i (V>W) = 0, 



(2.6) 



z— x,y 



which is linked to the unitarity property of the scattering 
matrix. Equation (J2.3I) then allows to describe the fol- 
lowing setups for the TI surface state. First, an orbital 
magnetic field has only effects when it is oriented perpen- 
dicular to the surface, B rb — B z (r,(f))e z . In cylindrical 
coordinates, we can then choose some gauge for the vec- 
tor potential A such that 



B z (r,, 



1 



{d r {rA^) - d^A r ) , 



(2.7) 



while A z drops out and is put to zero. Second, to describe 
the coupling of surface Dirac fermions to an in-plane ex- 
change field H(r) = (H Xl H y , 0), e.g., due to the magne- 
tization of a nearby FM layer, we write A = (c/cvf)M 
with M = (—H y ,H x ,0), where we used the spin Pauli 
matrices in Eq. (|2.5j) . For a Zeeman field, we can proceed 
in complete analogy where H now denotes the Zeeman 
field. A Zeeman or exchange field oriented along the e z 
direction can open a gap in the spectrum, and here we 
assume that such fields are not present. While the or- 
bital field breaks time reversal invariance, the Zeeman or 
exchange fields represent a time-reversal invariant per- 
turbation. Then A is determined by the magnetic field 
itself, and hence is not a gauge field anymore. 

Next we turn to graphene, where the Pauli matrices a 
are related to the two triangular sublattices constituting 
graphene's honeycomb lattice. We assume that no spin- 
flip mechanisms are relevant, i.e., physical spin is con- 
served. We can then focus on one specific Dirac fermion 
flavor with fixed valley index and spin direction. This 
excludes exchange or Zeeman fields, i.e., we put M = 
and hence A = A for graphene. Note that Zeeman 
fields in graphene are generally small compared to or- 
bital fields^ Moreover, we consider only smoothly vary- 
ing vector potentials such that it is indeed sufficient to re- 
tain only one K points Equation (12.3[) can then describe 
the following cases. First, we may have an orbital mag- 
netic field, precisely as for the TI case. Second, pseudo- 
magnetic fields generated by strain-induced forces ^SSr— 
or by various types of defects, e.g., dislocations^ also 
correspond to a vector potential, where time reversal in- 
variance implies that A has opposite sign at the two K 
points. A(r) can then be expressed explicitly in terms of 
the strain tensor^ where the resulting pseudo-magnetic 
field is also oriented along the e z axis and A z — 0. In 
addition, strain causes a scalar potential V(r), which is, 
however, strongly reduced by screening effects. The com- 
bination of orbital and pseudo-magnetic fields may al- 
low to design a valley filter, since the total (orbital plus 
pseudo-magnetic) fields can differ significantly at both K 
points^ 



B. Multipole expansion 

Our scattering theory approach considers magnetic 
fields [described by A(r) in Eq. (|2.3|) ] that smoothly 
vary on the scale of a lattice spacing and constitute a 
local perturbation, i.e., a well-defined cylindrical multi- 
pole expansion exists. Furthermore, A z = A, r = is as- 
sumed throughout. As we show below, for orbital fields 
we can choose a gauge where A r — 0. For strain-induced 
fields, strictly speaking, the problem is not gauge invari- 
ant, and we cannot impose gauge conditions. However, 
in a more narrow sense, a gauge degree of freedom still 
exists 4i For r — > oo, with complex- valued coefficients 

a im = ( a i-m ) > we then have the multipole expansion 



a$ 



**(r.*)-^+E £ 



im<$> 



J~ a l,m, 



(2.8) 



1=2 m=—oo 



where a denotes the total flux in units of the flux quan- 
tum $o = 27rfi,c/e. 

Let us now address the orbital magnetic field case, A = 
A, where we can exploit gauge invariance. We start from 
a more general situation with A r ^ 0, expressed as in 
Eq. (|2.8p with coefficients af^, and also allow for nonzero 
coefficients a^ m , . We now show that one can choose 

a gauge where A r = and a\ ' m , = 0. Indeed, gauge 
invariance implies that for arbitrary functions g(x, y), we 
are free to replace A; — > Ai + dig. Using a multipole 
expansion for rg(r, </>) with coefficients 5/ jTO , an equivalent 
gauge choice thus follows by the replacement 






Ar) 



imgi, m , 



We then choose the gauge function 



a 



(r) 



i« 



9l>i, 



l-l 



gi=i, m ^o 



In the new gauge, we arrive at Eq. (J2.8J1 plus the radial 
component 



A r = Y 



im4> 

\ r ) 



Using Eq. (|2.7p . the orbital field expansion (with r > 0) 
reads 



B z ( 



7\ = - 



oo oc 

E E • 

I— I rn— — oo 



rt+l 



(J-lJalS + trofc.xaii 



The to = term in A r neither generates flux nor mag- 
netic fields and can be omitted. Magnetic field profiles 
with af m 7^ arise only in time-dependent settings and 
will not be studied here. As a consequence, the radial 
component vanishes, A r = 0, and we arrive at Eq. 



III. SCATTERING THEORY 

For given energy E — hvpk, where k > through- 
out, the Dirac equation, Hip = Eip with Eq. (|2.3[) . has 
scattering solutions that we wish to obtain in the pres- 
ence of magnetic perturbations of the type in Eq. (J2.8I) . 
The solution for E = —hvpk follows simply by reversing 
the sign of the lower spinor component.— We are then 
looking for a solution ip{r, <p) — ipin + ?/Wt consisting, in 
the asymptotic regime r — > oo, of a plane wave (ex e tkx ) 
propagating along the positive x-direction, 



ipm(r,(f>) 



1 



a ikr cos (£> 



vr v 1 

plus the scattered outgoing spherical wave^ 



%ut{r,ct>) = F(<t>) 



e ikr j x 
oi<t> 



\f=2vr 



(3.1) 



(3.2) 



We adopt the same normalization conventions as 
NovikovJ^ From Eq. (I2.5[) we see that the incoming cur- 
rent density is j m = vf&x while (for the TI case) the spin 
density is (h/2)e y . Equation (|3.2p defines the scattering 
amplitude F((f>) for an outgoing wave deflected under the 
scattering angle <f>. The resulting scattered current den- 
sity implies the standard definitions of the differential 
(da/dtfi), total (otot) and transport (<7t r ) scattering cross 
sections ; 12 ' 58 respectively: 



S = TO)| a 



(3.3) 



o-tot = / # \F&)\ 2 = \/tt Im ^(0)- 
Jo V k 



<7tr = 



2tt 



(l-cos0) \F(4>)f 



The second equality for crtot expresses the 2D optical the- 
orem. When a random distribution of magnetic pertur- 
bations is present, the inverse mean free path determin- 
ing the conductivity is proportional to the transport cross 
section^ 



A. Born approximation 



that in real-space representation, the retarded Green's 
function (Hq — E 2 )" 1 is given by the Hankel function 



H 



(!) 12.60 
!~ 



i>out(r) = 



&J TlW 



d 2 r' H^> (k\r - r'\){cr ■ p' + hk) 



2V2h^ Q 

The asymptotic large-p behavior of the Hankel function 
(where 77 = 1,2 = ±) ia^ 



H<p{p)~J±e 

V ftp 



±i(p-(2i/+l)7r/4) 



(3.5) 



which implies that tjj ou t for r — > 00 indeed has the form in 
Eq. (|3.2p . After some algebra, we obtain the scattering 
amplitude in Born approximation, 



F{4>) 



/2irk 
$0 



l *' 2 / rdr / d(f>' 



S11K 



x e -2ifcr|sm(0/2)|sin0' A (r, </>' + 0/2). (3.6) 



For radially symmetric perturbations, d^K.^ = 0, the 0'- 
intcgration can be done, and we obtain 



F(0) 



-2m- 



Jl-Kk 



-ic/>/2 



(3.7) 



x / rdr Ji (2fcr|sin(0/2)|) A^(r), 
Jo 



with the Ji Bessel function. 



B. Radially symmetric case 

Next, we address the full (beyond Born approximation) 
scattering solution for radially symmetric perturbations, 
A = A^(r)e0. In that case, the total angular momentum 
operator J z in Eq. (J2.4I) is conserved and has eigenval- 
ues Hj with j = m + 1/2 (to 6 Z). We thus expand 
the spinor wavefunction in terms of angular momentum 
partial waves ip m (r) = {f m ,ig m ) T , 



For small perturbation A^r, 0), one can evalu- 
ate the scattering amplitude within the first Born 
approximation.— Strictly speaking, the long-ranged part 



A^ ex a/r in Eq. (|2.8[) can not be treated perturbatively, 
and in this section we assume a = 0. 

The unperturbed state is the incoming plane wave ip[ n , 
Eq. p. II) . Within lowest-order perturbation theory, the 
scattered wave obeys 



[H - E]ip out 



evF 



A,/ 



<T Ipi, 



(3.4) 



where Hq is the unperturbed Dirac Hamiltonian. Mul- 
tiplying both sides of Eq. (|3.4[) by H + E and noting 



ip(r, 0) 



_ p -i<j»J z /2 V^ j(m+l/2) 



m— — cxd 



where the Dirac equation yields 



• \ d r + —\ cf x 



m + 1/2 + (p(r) 



ipm(r), (3- 



1p m = klp n 



The magnetic flux (in units of the flux quantum $o) en- 
closed by a circle of radius r around the origin is 



tp(r) 



2nr 



Mr), 



(3.10) 



where a = ip(oo) in Eq. (|2.8[) . The continuity relation 
(|2.6[) must hold for each partial wave tp m separately, and 
implies 



d r (r^ m a x tp m ) = 0. 



(3.11) 



Introducing dimcnsionlcss radial coordinates, p = kr, a 
closed equation for the upper component, f m (p), follows, 



- p d p ( P d p ) + i-(J^ + w m + w^ 



W m (p) = 



1/2 + <p(p/k) 



fm = 0, 

(3.12) 



where W m = d p W m . The lower component is obtained 
from 



9m{p) 



d P + Y~ W ™) /< 



(3.13) 



These relations imply a general expression for the scatter- 
ing amplitude F(<f>) under radially symmetric magnetic 
perturbations, and thus for the various cross sections in 
Eq. (|3.3|) . For p —> oo, the term ex a/r in Eq. (|2.8[) dom- 
inates and the general solution to Eq. (|3.12[) is given in 
terms of Hankel functions, 

fm{p) = a m H^ +a (p) + b m HW(p), (3.14) 



As a consequence, qualitatively different effects beyond 
the electrostatic case arise, such as the AB effect. The 
cross sections in Eq. (|3.3j) are then given by 



CTtot 



<7 tr 



k^ 

m 

k^ 



sin 2 (S m ) . 



sin 2 (S m+ i - 6 m ) . 



(3.19) 



Scattering theory has thus been reduced to the determi- 
nation of the phase shifts S m . In the electrostatic case^ 
the phase shifts obey the symmetry relation S m = 5- m -\, 
implying the absence of backscattering, F(ir) = 0. In the 
magnetic case under consideration here, in general this 
symmetry relation breaks down, and hence backscatter- 
ing is not suppressed anymore, F(ir) ^ 0. This is closely 
related to the fact that magnetic fields can confine mass- 
less Dirac particles. 32 



IV. MAGNETIC DIPOLES 

As a first application, we analyze the scattering of 2D 
massless Dirac fermions by a fixed magnetic dipole mo- 
ment m located at position r = (0, 0, h), i.e., at a height 
h above the origin of the 2D plane. In that case no total 
flux is generated, a = 0. The results of this section are 
obtained under the Born approximation, see Sec. MI Al 



with complex coefficients a m and b m . The lower spinor 
component then follows from Eq. p. 131) . 

9m(p) = a m H^ } +a+1 (p) + b m H^l a+1 (p). (3.15) 



2i<5„ 



The continuity relation (|3.1ip implies a m = b m e z " Jm , i.e., 
the outgoing wave can differ from a free spherical wave 
only by a phase shift S m , which depends on the magnetic 
perturbation and is determined in Sec. |Vj Using the 
Bessel function expansion formula 



%p COS 



E 

mSZ 



■m imd> 
I e v 



Jm(p) 



r(1.2) 



and the asymptotic behavior of Hi ' , see Eq. (|3.5j) . we 
find 



L _ -ina/2 

2 



(3.16) 



We then obtain the scattering amplitude in terms of 
phase shifts as for the electrostatic case^ 



F(4>) 



/Wk 



J2 (e 2l5m ~ 1) e im *, 



mGZ 



but 5 m includes the total flux a, 



ira/2. 



(3.17) 



(3.18) 



A. Perpendicular orientation 

First, we consider a dipole moment oriented perpen- 
dicular to the layer, m± = m±e z , where we have the 
isotropic (vector potential) perturbation 



A^(r) = m± 



(r 2 + ^2)3/2 



(4.1) 



The Born approximation, see Eq. Q3.7p . yields the scat- 
tering amplitude 



F ± ((f>) = -ie-^ /2 y/{2ir) 3 k (m J 7$ o )e~ 2 ' c,l|sin(0/2)l , 

(4.2) 
and the transport cross section is 

<Ttr,a. = (2^) 4 fc (m ± /$ ) 2 F 2 (fc/ 1 ). (4.3) 

Here we define the functions 

j.n a —4:tx 



4 f 1 t Tu e 

F n (x) = - dt - 7 ==: 

7T Jo Vl - t 2 



(4.4) 



which can be expressed in terms of hypergeometric func- 
tions. As shown in Fig. [TJ after reaching a maximum 
around kh ~ 0.31, the transport cross section (j4.3[) de- 
creases with increasing energy and approaches zero for 
E — » oo. Figure [T] also shows a polar graph for the dif- 
ferential cross section, da/d(f> — |F(0)| 2 , at kh = 0.2. 
Evidently, scattering by a dipole oriented along the e z 
axis is almost isotropic. 




FIG. 1: (Color online) Born approximation results for the 
transport cross section a tr for scattering of massless Dirac 
fermions on a magnetic dipole. We show ctr in units of 
(27r) 4 (|m|/$o) 2 vs energy E in units of Hvf/Ii, where h is 
the distance of the dipole from the layer. The dot-dashed 
blue curve is for the perpendicular orientation, see Eq. (|4.3p . 
the other curves are for the parallel orientation and several 
angles rj (cf. legend), see Eq. H4.8[) . Upper right part: Po- 
lar plot of the differential cross section, |_F((/>)| 2 , for E = 0.2. 
Different curves correspond to the ones in the main panel. 



C. Bilayer graphene 

Let us briefly comment on the results under a quadratic 
dispersion relation as realized in bilayer graphene^ Re- 
peating the Born approximation analysis, the scattering 
amplitude is found to contain an additional cos(</>/2) fac- 
tor modifying the above expressions. In fact, we find 
instead of Eqs. fj4T3|) and (|H5]i : 

a[f L ± G) = (27r) 4 fc (m ± /$ ) 2 (F a - F 4 ), 

a^ G} = (27r) 4 fc (mn/$ ) 2 (4.9) 

x F 4 - F 6 + (F a - 3F 4 + 2F 6 ) cos 2 (rj)] , 

with F n — F n (kh). The quoted expressions hold for the 
quadratic dispersion relation of bilayer graphene, and co- 
incide with the results for the conventional Schrddinger 
case. In contrast to the monolayer results for cr tr ii in 
Eq. (|4.8j) . the transport cross section (|4.9j) carries no 77- 
dependence at low energies. The latter is a distinctive 
feature of 2D massless Dirac fermions. Finally, we note 
that the cr t r,|| results for 77 — in Fig. Q] coincide with the 

bilayer result <7 tl . j_ (when m± = mi) for perpendicular 
orientation. 



B. Parallel orientation 



RING-SHAPED MAGNETIC FIELDS 



If the dipole instead points parallel to the layer, 

m = m\\ [— sin(?7)ea; + cos(rj)e y ] , (4.5) 

where 77 denotes an angle, we have 



In this section we consider the scattering states for 
a radially symmetric ring-shaped magnetic field. The 
scattering setup is schematically sketched in the inset of 
Fig. El 



A (r,0) 



II 



h M4>-v)y {r2 + h2)3/2 r 



, (4.6) 



i.e., no radial symmetry is present. Using Eq. (|3.6|l . we 
now find the scattering amplitude 



*||(0) = -e-^/ 2 N /(27J3fc(m||/$o) (4.7) 

x cos(77-(/./2)e- 2 ' £ ' l|sin(</,/2)l , 

and thus the transport cross section 

a tl ,H = (27r) 4 fc (m||/$ ) 2 [f 4 + (F 2 - 2F 4 ) cos 2 (77) 

(4.8) 
where F„ = F n (kh). This cross section depends on 
the orientation rj of the dipole even for small energies, 
kh <C 1. When averaging over rj (which is equivalent to 
setting 77 = vr/4), we find er tr ji = (mn/m±) 2 <7ti-,±/2- The 
transport cross section (14.81) has a maximum for an in- 
dependent energy, see Fig. [T] and again approaches zero 
for E — > 00. The differential cross section shown in Fig.Q] 
also reveals a pronounced angular dependence tied to the 
orientation of the dipole. 



A. Innnitesimally thin ring 

Let us first study the exactly solvable model of an 
innnitesimally thin ring of radius R around the origin, 
where A^r) follows from Eq. (|3.10[) with 



ip(r) = aO(r — R), 



(5.1) 



where O is the Heaviside step function and, as before, a 
is the dimensionless total flux through the ring surface 
area. For the orbital field case, this implies B z (r) — 
(a^ /2nR)S(r - R). With p = kr and K = kR, the 
solution to Eq. (f3~T2"l) is 



fm(p) = 



(e 21 ' 



iJm(p), 



5 m rj-(l) ( n \ 
H m+a(P) 



H-m+a\P) 



p<TZ, 
P>K, 

,(5.2) 

with b m in Eq. (|3.16|) . The requirement of continuity of 
tp m (r) at r = R, together with Eq. (|3.13p . leads to two 
boundary conditions for f m . With 7?. ± = 1Z ± + , they 
read 

f m (n + ) = f m (K-) , f m (n + )-f m (n-) = ^f m (n), 



n J 



(5.3) 




FIG. 2: (Color online) Transport cross section a tI (in units 
of 2/fc) vs dimensionless flux a for a finite-width magnetic 
ring, see Sec. IVCl with kRi = 0.01 and R 2 = 2i?i. (We 
here also allow for a < 0.) The numerical results are close 
to the ideal AB prediction for the infinitely thin solenoid, 
(Jtr = (2/fe) sin 2 (7ra). Inset: Schematic scattering geometry. 
The plane wave (blue solid arrows) coming in along the e x 
direction is scattered by a ring-shaped magnetic field present 
for Ri < r < R.2 (shaded region). The outgoing spherical 
wave is indicated by red dashed arrows. 



where again /' = d p f . The coefficient a m and the phase 
shift 8 m appearing in Eq. (|5.2p then follow from the 
boundary conditions (|5.3[) . In particular, when J m (lZ) ^ 
0, the phase shift S m can be determined by evaluation of 
the logarithmic derivative 

^ ^ d]nf m (p = K+) = a + J' m {K) 



dp U J m {TL) 

m + a J m+ i(TZ) 



n 



J m (TZ) 



where we used the second boundary condition in 
Eq. (|5.3j) . As a result, with the Neumann function Y v , 
we find 



tan(5„ 






while a m is given by 



r(2) 



i = b 



J m (n) 



(5.5) 



(5.6) 



Equation (|5. 5[) stays valid beyond the thin-ring limit 
when a more general form for C m is used, see Sec. IVCl 
For the special case J m (7£) = 0, Eq. (|5.2|) implies 
e «? m = -tf$ (ft)/fl$. a (ft) and, using /£,(*+) = 



r(l) 



(2) 



am ~ m JUn) 

Equations (J5.5I) and (15.61) include these relations when 
taking the limit J m (lZ) — > and C m — > oo. 



B. Aharonov-Bohm scattering amplitude 

Let us first consider the R — > limit of the above set- 
ting, which corresponds to the pure solenoid case. This 
allows us to study the Aharanov-Bohm (AB) effect for 
ultra-relativistic Dirac fermions. In order to extract the 
singular part, we first rewrite / TO (r) in Eq. (13. 14)) as 



f m (r) = 2b m — 



sin(7ra - 5 m ) J m+a (kr) 



sin(7ra) 
+ (-) m sin(<5 m ) J_ (m+Q) (fcr) 



(5.7) 



Imposing regularity for f m (r) as r — > requires the phase 
shift (|3.18[) to be S m — — (7ra/2)sgn(m + a). Correspond- 
ingly, for i? — t> 0, the scattering amplitude p,17[) is given 
by 



F(4>) = 



(e-*™-!) £ 



o ™0 



(5.8) 



W-i 



ra _ l} ^ 



im^i 



?n— — oo 



where a = [a]+{a}, with integer part [a] and non-integer 
part < {a} < 1. Summation of the series in Eq. (|5.8|l 
yields^ 



F(0) = 



:(2TrS(<j)){cos(TTa)- 1] (5.9) 



+ r -i(M+i/2)0 sin(Tra) \ 
sin(<£/2)/' 



Up to the forward scattering (<f) = 0) amplitude, Eq. (|5.9[) 
reproduces the AB result t 45 i 46 i 48 here obtained in terms 
of scattering phase shifts. Note that the forward scat- 
tering (5-term, missing in the AB calculation^ naturally 
appears in our phase shift analysis and is essential for 
establishing unitarity of the scattering matrixi 47 ' 52 

In alternative approaches to obtain F((j>) for the ideal 
solenoid, following the original AB method^ the asymp- 
totics of the exact wavefunction is computed from its 
integral representation. As a result, the incident wave 
corresponding to Eq. (|3. 1 1) has an additional phase fac- 



tor e~ ma sgn(sm0) e -;a^ ie ^ one has a mu i t i- va lued in- 
coming plane wave. The precise relation between these 
two approaches has been discussed in several works and 
is still under debate^Zr— albeit the difference is of lit- 
tle relevance to experimentally observable quantities. In 
particular, the transport cross section er tr in Eq. (|3.3[) 
does not depend on the forward scattering amplitude at 
all. We conclude that our approach is able to reproduce 
the AB effect, <7t r = (2/fc)sin 2 (-7ra), with oscillations as 
function of the dimensionless flux parameter a. In par- 
ticular, <7 t r = for integer a. 



C. Magnetic ring of finite width 

Before discussing concrete results for the scattering 
amplitude and the transport cross section in the pres- 
ence of a ring-shaped magnetic field, we now generalize 
the setup to a finite width, with R\ < R 2 denoting the 
inner and outer radii of the ring, cf. the inset of Fig. [2] 
Again, A^(r) in Eq. (|3.10[) is expressed in terms of a 
dimensionless flux function ip(r). When A^ is a vector 
potential, the associated magnetic field B z (r) = B is 
taken uniform within the ring region and zero outside; 
for concreteness, we take B > 0. This profile allows for 
an exact solution, while more general smooth field pro- 
files can be treated within the Wentzel-Kramers-Brillouin 
(WKB) approximation, see Sec. IV Dl 

We use dimensionless coordinates (p = kr and 72.1.2 = 
kRi 2 ) and flux parameters, 



V\,2 



nBR 2 12 



^2 = :^2 > a = v 2 -v x . (5.10) 



The function <p then reads with r = p/k: 



<p(r) 



0, 
vp 2 - u u 

a, 



P<n ll 

fti <p<n 2 , 

p>u 2 . 



(5.11) 



For i?i — > R 2 , this reduces to Eq. (|5.1[) . In particular, a 
in Eq. (|5.11[) again denotes the total dimensionless flux. 
For given j — m + 1/2, the components of the Dirac 
spinor ip m obey Eqs. (p7l2|) and (|5TT5)| . with Eq. ([STTTj) 
now determining W m (p). The solutions for r < R% and 
r > R 2 are as in Eq. (|5.2I) , 



fm{p) 



r(D 



iJm(p), p<1Zi, 

[e 2l6m H^+c (p) + Hi% a (p)), p>K 2 , 

(5.12) 
where a m and S m are to be determined, and b m is given 
in Eq. (HTTp) . For i?i < r < R 2 , Eq. ^53g) can be solved 
in terms of the confluent hypergeometric functions <£> and 



fm{p) 



6, 



p \™\ e -vp 2 /2 



Cm*(Cm,l+ |m|;^p 2 ) 



+ d m *(Cm,l + |m|;^ 2 ) 



(5.13) 



1 + rh<d(m) — 1/4^, m = m — v\. 



The coefficients c m and d m , together with a m and the 
phase shift 8 m in Eq. (|5.12p . follow by matching ip m at 
r = Ri=i j2 . Taking into account that W m is a continuous 
function of p, we have 



f m hit) = U (n; 



fin fit) = fin (^) . ( 5 - 14 ) 



where the second condition follows by continuity of the 
lower spinor component g m . It is convenient to introduce 
the transfer matrix T m connecting the solutions at p = 








FIG. 3: (Color online) The main panel is as in Fig. [2] but for 
kRi = 1.32 and R2 — 5-Ri. Inset: Current density j${p) = 
e<t> ' 3 vs radial coordinate p for the quasi-bound state with 
j = 3/2 present at a ~ 6. 



72|~ and 72 



' m 



*^m I m 



2 i 

Jm(72l) 

■4(^0 

(5-15) 
Explicitly, the transfer matrix for the magnetic ring of 
finite width is 



Jm ( 'v-2 
fin(^2 



f m (nf 

fin(R>t 



I m 



$ 2 $ 2 



*a *2 A *i *i 



$1 \I>i 



(5.16) 



where we use the abbreviation 

$ i= i, 2 = TZ^e-^HiU, 1 + WW), 

and similarly for ^. We mention in passing that for 
the infinitesimally thin magnetic ring in Sec. IV Al (where 

i?i = R 2 = R), the transfer matrix is i m 



a/72 1 

For the finite-width ring, using Eq. (|5.12[) the phase 
shift 5 m is then again given by Eq. (|5.5p . with 72 — > 722 
and the logarithmic derivative C m replaced by 



Mm, 2 



«m,l 
»m,2 



— /m 



Jm (72l) 
Jm(72l) 



(5.17) 



With the above expressions, it is straightforward to 
compute the scattering phases 5 m = 8 m — -kol/2 numer- 
ically for the finite-width ring geometry. Thereby we 
obtain^ the scattering amplitude F{<f) from Eq. (13.17)) 
and the transport cross section <7t r from Eq. (|3.19[) . 

Numerical results obtained under this approach are 
shown in Figs. [2] and [3] First, in the main panel of Fig. [51 
we show the transport cross section er t r as a function of 
the total flux a. In this example, both radii R\ and 
R 2 were chosen very small, such that scattering by the 
ring is close to the one by an ideal AB solenoid. As a 
consequence, we observe the AB oscillations with unit 



flux period. In contrast to the ideal AB result, a com- 
plete suppression of scattering for a G Z is observed in 
the finite- width ring only for a = 0, while the maximum 
value CTt r = 2/k for half-integer a is still perfectly real- 
ized. In fact, the phase shift analysis in Sec. IV Bl shows 
that a given oscillation period is determined by one spe- 
cific m value in the ideal AB case. For the non-ideal 
finite- width ring, other total angular momenta also start 
to contribute, and this mixing effect destroy the perfect 
constructive interference needed for tr tr = 0. On the 
other hand, the destructive interference responsible for 
the maxima of trtr at half-integer a is more robust since 
it is dominated by a single m value. 

In Fig. [3l we study scattering by a much larger ring. 
In this case, the AB effect is absent, which can be under- 
stood by noting that the Fermi wavelength (2%/k) of the 
particle is now smaller than the outer circumference lixRi 
of the ring. Quantum interference of waves surrounding 
the obstacle in opposite directions is then largely aver- 
aged out, and, moreover, the wavefunction can partially 
penetrate into the ring area. However, a remarkable peak 
feature at a w 6 appears now in the transport cross sec- 
tion, see Fig. [3] This feature can be traced to the appear- 
ance of a quasi-bound state with j = 3/2 at this flux value 
(for the considered energy), which then causes a scatter- 
ing resonance, cf. our discussion in Sec. lVDl The inset of 
Fig. [3] shows the current density profile for precisely this 
quasi-bound state. While the radial component vanishes, 
j r = 0, we find a circularly oriented current, j^ =/= 0, 
which is mainly localized inside the ring (r < R\) and 
represents a current-carrying bound state. We note that 
more quasi-bound states appear for larger a, causing ad- 
ditional peak features in <J tr (a) beyond those shown in 
Fig.H 



D. Quasi-bound states and scattering resonances 

The magnetic confinement built up by the ring-shaped 
field can generate quasi-bound states, which for R2 — > 00 
become true bound statesj 32 i 33 The quasi-bound state 
spectrum then causes resonances in the scattering am- 
plitude when the energy E = hvpk is varied. For given 
total angular momentum j = m + 1/2, the correspond- 
ing phase shift S m (E) goes through the value n/2 as E 
crosses a resonance level E r . The corresponding reso- 
nance width r r can be estimated frorn^ 

(d/dE) cot [5 m (E = E r )] = -2/T r . 



To access these resonances, we first put Eq. (|3.12|) into a 
canonical form with separated kinetic and potential en- 
ergy terms. The substitution f m (p) = P 1 ^ 2 fm{p) yields 

[-d 2 p + v m ( P )] f m = ; m , v m = w 2 m + w; n , (5.18) 

where V m (p) is an effective potential energy for the ra- 
dial motion and the lower spinor component is g m = 
p' x / 2 {-d p + W m )f m . In this form, Eq. (1535)) can be 
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FIG. 4: (Color online) Partial cross section sin 2 5 m vs energy 
E for a ring-shaped confinement as in Sec. lVCl The numerical 
results are for total angular momentum states with m — 1 
(solid black) and m — —1 (dashed blue curve). The radii 
are R2 = 7Ri and Ri — 0.5£b with £b — \/2c/eB, and 
energies are in units of Kvf /&b- Inset: WKB results for quasi- 
bound state energies E r vs R2 (lengths in units of £b), for 
m = 1 (black circles) and m = — 1 (blue diamonds) with fixed 
Ri = 0.51b- For comparison, the exact levels for infinite R2 
from Ref. [33 are shown for m — 1 (dotted black) and m — — 1 
(dashed blue curve). 



treated within the standard WKB approach, which repre- 
sents an attractive alternative to semiclassical approaches 
to the Dirac equation as it avoids the appearance of non- 
Abelian Berry phases! 42 ' 62 ' 63 For a magnetic ring as in 
Sec. lVCl the effective potential V m has a hard repulsive 
core for r — > plus a barrier at larger distances, i.e., a 
quantum well is formed with classically allowed motion 
for ro < r < r\. The "turning points" ro,i here depend 
on the energy E — hvpk under consideration. For fi- 
nite i?2, this barrier is of finite width and quasi-bound 
states within the well region may exist. The classically 
forbidden region r± < r < r^ (where ri is another turning 
point) then corresponds to tunneling trajectories where 
the "particle" escapes from the well region. For R2 — » 00, 
the barrier becomes infinitely wide and this escape prob- 
ability vanishes, i.e., we obtain true bound states in the 
well region. Using the radial variable r = p/k, Eq. (|5.18[) 
reads 



[-d 2 + U m (r)]f m (r) = ef m (r), 



(5.19) 



U m (r) 



u> 



+ d r w m , w m (r) = kW m (kr), 



10 



where the modified Bohr-Sommerfeld quantization con- 
dition for the complex-valued "energy" e = k 2 i 



,64 



dr y/e- U m (r) = n ( n + - - 

r \ Z 

i / r(io + i/2) 

x(o)= « In U-*o + l/2)[l + « 

+ a(l — In a), 
1 m 



XJa) 
2ir 



,(5.20) 



-27?a] 



• i = - I dr \/U m (r) - e, 

7T 



with n = 0, 1, 2, . . . and the Gamma function T(z). The 
complex resonance values for e solving Eq. (|5.20[) can 
be found numerically. Equation (|5.20[) is formally exact 
for the case of a parabolic barrier, but also applies for 
an arbitrary smooth potential and is expected to remain 
accurate" even for small n. We now write 

e = (k — ij/2) sa k — ikj. 

For a quasi-bound level with energy E r — hvpk, the reso- 
nance width is then T r = hvpj. Using Im%(a) ss e~ 27ra /2 
for a > 1, we obtain 



T r /h- 



TfV 



2-na 



(5.21) 



where the period of radial motion is 



2fc 

T k = — 



dr 



vf J ro y/k 2 - U m (r) 

Our numerical results for the partial cross section, 
sin 2 5 m , as a function of energy, and the WKB results 
for the corresponding quasi-bound state energies E r are 
shown in Fig. 0] Here we take the field profile as in 
Sec. IV Cl — With increasing i?2 (keeping i?i fixed), new 
quasi-bound energy levels localized in the well region ap- 
pear, see inset of Fig. |4] Very good agreement with ex- 
act quantum calculations 33 for the infinite barrier case 
(i?2 — > oo ) is observed, i.e., these energy levels remain ba- 
sically unchanged when increasing i?2- The only notice- 
able deviation from the exact spectrum of Ref. l32| is seen 
for to = 0, where the potential U m —o(r) creates an in- 
finitely attractive well for r — > 0. In that case, the WKB 
approximation becomes questionable in that "steep" re- 
gion. The main panel in Fig. U illustrates the sequence 
of quasi-bound states present because of the magnetic 
confinement. The corresponding scattering resonances 
appear as peaks in the transport cross section ctr when 
varying energy or the effective flux parameter a. 



VI. CONCLUDING REMARKS 

In this paper, we have studied scattering of massless 
two-dimensional Dirac fermions by magnetic perturba- 
tions of various types. The model is applicable to quan- 
tum transport in monolayer graphene and for the sur- 
face state of strong topological insulators. The magnetic 
fields can correspond to orbital or Zeeman fields, strain- 
induced fields in graphene, or exchange fields generated 
by ferromagnets. 

The full scattering solution was discussed in detail for 
radially symmetric perturbations, where the scattering 
amplitude can be expressed in terms of phase shifts in a 
given total angular momentum channel, and within the 
Born approximation for the general case. Our approach 
now allows for a systematic study of the scattering of 
Dirac fermions on magnetostatic perturbations. 

As applications, we have studied scattering by mag- 
netic dipoles within the Born approximation, and fully 
nonperturbative scattering for the case of ring-shaped 
magnetic fields. The Born approximation is only valid 
when the perturbation has zero total flux (a — 0). For 
the magnetic dipole, we have pointed out characteris- 
tic angular dependencies in the differential cross section 
that may allow to unambiguously identify massless Dirac 
fermions. For the ring-shaped field case, as one increases 
the lateral size (R2) of the magnetic perturbation, we 
have a crossover from the Aharonov-Bohm case to a 
regime dominated by scattering resonances. In the first 
case, R2 — > 0, particle trajectories surround the flux re- 
gion but essentially do not penetrate it, leading to the 
oscillatory transport cross section er t r oc sin (na). In 
the second case, where the particle wavelength is small 
against the size of the perturbation, fci?2 > 1, the AB 
oscillations in Ot r (a) are absent. However, now quasi- 
bound states arise due to the magnetic confinement, 
causing scattering resonances which show up as peaks 
in <7 tr (a). 

To conclude, we hope that these predictions motivate 
further theoretical work and that they will be tested ex- 
perimentally in the near future. 
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